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We study gravitational waves to first and second order in amplitude in vacuum asymptotically flat 
spacetimes. The Einstein equations are solved to first order and these solutions are superposed to 
form a time-symmetric ingoing and then outgoing pulse regular everywhere. The waves are assumed 
to have odd-parity and a non-vanishing angular momentum which keeps them away from the axis 
at all times. The averaged energy of the waves is evaluated. The relevant Einstein equation is then 
solved to second order in the amplitude. The influence of the angular momentum of the waves on 
the rotation of local inertial frames with respect to the frames at great distances is analyzed. The 
rotation of the frames occurs even in the region around the origin where spacetime is almost flat. 
The rotation is without time delay as it follows from the constraint equation. 

The results are illustrated graphically for various values of the "harmonic index m" corresponding 
to azimuthal rotation and the "harmonic index V describing the latitudinal rotation of the waves. 
The apparent motions of the fixed stars on the celestial sphere as seen through rotating waves from 
the local inertial frame are calculated and displayed. 
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I. INTRODUCTION 

It was just 100 years ago in Prague when Einstein 
wrote the paper [l[ in which he, for the first time, ex- 
pressed his understanding of Mach's Principle. Within 
his pre-general relativity theory in which there was only 
one metric function he considered a mass point inside 
a shell accelerated "upwards" and found that the mass 
point is dragged along by the shell. 

Many formulations and studies of Mach's Principle 
have appeared during the last 100 years, most of them 
were analyzed in the Tiibingen conference in 1993 which 
led to the remarkable volume @ containing lectures as 
well as valuable discussions. We studied Machian effects 
in various contexts, both in asymptotically flat space- 
times and within cosmological perturbation theory - see, 
e.g., Q, and number of references therein; later, cf. 
Schmid 0. 

More recently, we investigated the subtle question 
whether dragging of inertial frames should also be at- 
tributed also to gravitational waves. After the discovery 
of binary pulsars losing energy and angular momentum 
as a consequence of emitting gravitational radiation it 
would be surprising if gravitational waves did not have an 
influence on local inertial frames. However, there are still 
doubts about the status of gravitational stress-energy as 
compared with stress-energy tensor of matter in re- 
lation to Machian ideas (see, e.g., jlj, p. 83). 
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In [5| and [6(, we studied gravitational waves with 
translational symmetry that rotate about the axis of 
cylindrical polar coordinates. We solved the Einstein 
equations to first order in the wave amplitude and formed 
a pulse. By solving the relevant Einstein equation to sec- 
ond order we calculated and illustrated graphically the 
rotation of a local inertial frame at the axis which is 
caused by the rotating wave pulse. 

In the present work we investigate the effects of rotat- 
ing gravitational waves in a more general, asymptotically 
flat setting, without assuming any symmetry. We again 
start out from linearized theory and construct an ingo- 
ing rotating pulse of radiation which later transforms into 
an outgoing pulse. While in the translational symmetric 
case our waves were characterized by just one harmonic 
index m governing the number of wave crests in tp, now 
the situation becomes considerably richer involving both 
spherical harmonic indices I and m. 

Our problem now resembles, in some respects, the 
study by Lindblom and Brill Q of the gravitational 
collapse of a slowly rotating, massive spherical shell 
freely falling under its own gravity. We reconsidered 
and extended their analysis in Q. The construction 
of our waves starts out from (the linearized version of) 
a Bonnor- Weber- Wheeler time-symmetric, ingoing and 
outgoing pulse of a typical width a and a typical am- 
plitude C. One can make a pulse that is rotating, but 
still localized in the radial direction so that it resembles 
a falling and rotating shell. Although we do not find 
the complete metric up to the second order in the am- 
plitude, the resulting spacetime will be asymptotically 
flat at spatial infinity since the first-order perturbations 
decay as h\}„ ~ r _ ( i+1 ) there, / being the multipolar in- 
dex. Hence, one can define "fixed stars" at infinity at rest 
in an asymptotically Minkowski system given in asymp- 
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totically spherical coordinates t, r, 9, ip. The waves have 
mass-energy, a nonvanishing total angular momentum, 
and a vanishing total linear momentum. Since we con- 
sider weak waves it is plausible to assume that a global 
coordinate system can be extended over the whole space- 
time. At infinity it becomes the rest frame of the fixed 
stars. 

Near the origin, the first-order metric of our waves be- 
haves as r l , so the region around the origin will be very 
nearly flat for I sufficiently large. When, however, a local 
inertial frame is introduced at the origin, we find that 
its axes rotate with respect to the lines ip = const of the 
global frame, i.e., with respect to fixed stars at infinity. 

Near the origin the congruence of the worldlines ip — 
const twists and the observers attached to these lines ex- 
perience Euler acceleration proportional to d t uj , where 
ujq is the angular velocity of the inertial frame near 
the origin. The angular velocity luq enters the second- 
order odd-parity dipole I = 1 perturbation of the met- 
ric, g$ = —ujq r 2 sin 2 9. (The centrifugal acceleration 
is higher order in the angular velocity.) The situation 
thus indeed resembles the interior of a collapsing slowly 
rotating shell - see Q where the vorticity of the lines 
ip = const is given in covariant form. In [8[ we also cal- 
culated how the fixed stars at infinity rotate with respect 
to the inertial frame at the origin by considering photons 
emitted radially inwards from the stars. 

In the present work we first consider gravitational 
waves with odd parity in the Regge- Wheeler gauge in 
the linearized Einstein theory. The metric components 
h^u describing the waves can all be derived from one 
potential-type function satisfying the flat-space scalar 
wave equation. We choose this function in the form of a 
generalized Bonnor- Weber- Wheeler ingoing and outgoing 
pulse characterized by specific harmonic multipolar or- 
ders I and to, | to | < I. For m ^ 0, the field represents the 
superposition of rotating waves which are ingoing from 
infinity towards the origin, stay regular all the time, and 
become outgoing waves again to decay at infinity. 

(2) 

Turning to the second-order metric h^J , we employ the 
second-order perturbative Einstein equations which are 
well known to be again linear differential equations that 
contain as "source terms" expressions quadratic in the 
first-order perturbations h$ and their first- and second- 
order derivatives (see, e.g., [HI)- Since we look for 
the second-order rotational perturbations we concentrate 
on odd-parity again. To determine the influence of grav- 
itational waves on the rotation of local inertial frames at 
the axis of symmetry we do not need to solve the second- 
order Einstein equations in general, because the high I 
terms do not have significant amplitude at the origin. 
The terms that give the inertial frame near the origin 
are I = l,m = 0, so we can average the "source terms" 
over azimuth ip, and the problem becomes axisymmetric. 
This simplifies the calculations. In addition, an impor- 
tant result proved and discussed in [ll|, [H[ states that 
the Regge- Wheeler gauge can be introduced in both first 



and second order, and it is unique. As emphasized in 
[lH, perturbations in the Regge- Wheeler gauge can be 
considered to be gauge-invariant expressions. 

The angular velocity wo of rotation of the local iner- 
tial frames at the origin, whose neighborhood remains 
as flat as we wish for all times by choosing very high I, 
is found by solving a simple linear elliptic equation for 
the second-order dipole odd-parity perturbation with a 
complicated source term given by the averaged compo- 
nent, (Rty [h^\ h^}), of the Ricci tensor, quadratic in 
h^ 1 ' and its derivatives. 

A lengthy, complicated form of (RlJ) allows us to find 
wo explicitly in a useful form only at particular values 
of time. However, by using computer algebra system we 
produce nice "bell-shaped" curves showing how uj in- 
creases as the waves approach the origin and decreases 
as the waves recede back to infinity. From loq plotted 
for different azimuthal harmonic index m, with I fixed, 
it is easily seen how an increase of the azimuthal angu- 
lar momentum of the waves implies a stronger dragging 
of inertial frames (gyroscopes) in nearly flat regions sur- 
rounded by the waves. 

We also calculate how the celestial sphere will look 
when observed from the inertial frame at the origin. 
This is influenced both by the passage of the stellar light 
through gravitational waves and by the dragging of the 
inertial frame due to their rotation. The resulting figures 
describing the apparent motion of fixed stars on the ce- 
lestial sphere nicely illustrate the role of the multipolar 
harmonic indices I and to on the pattern of the motion. 

The paper is organized as follows. Since odd-parity 
perturbations follow from a function satisfying the wave 
equation, in the next Sections II and III we study rotat- 
ing scalar-field waves in Minkowski space. We evaluate 
their energy density and angular momentum density av- 
eraged over the azimuth ip and in every direction, after 
averaging also over the latitude 9. We turn to rotating 
linearized gravitational waves in Sections IV and V. Re- 
calling some basic properties of the Regge- Wheeler for- 
malism and tensor harmonics, and fixing our convention, 
we calculate the averaged energy density of the rotating 
linearized gravitational waves and compare it with that 
of scalar waves. We find that for "rapidly rotating" waves 
with to ~ I the scalar field part dominates. 

The core of the paper lies in Sections VI and VII. Here 
we consider the second-order perturbations. After ex- 
panding in tensor spherical harmonics, and taking the 
axially symmetric part, we find the second-order Einstein 
equations for the radial parts of the second-order odd- 
parity metric components in the Regge- Wheeler gauge. 
For general I we get linear wave-type equations with 
source terms given by combinations of averaged compo- 
nents (R^u [h^, h^]} . For the dipole perturbations we 
obtain an elliptic equation which we solve by variation of 
constants. The angular velocity loq of an inertial frame 
near the origin is determined and analyzed both ana- 
lytically and numerically in Section VII. Here also the 
figures demonstrating the dragging of inertial frames by 
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the waves and the apparent motion of fixed stars as seen 
from the inertial frame at the origin are presented. In the 
concluding remarks we discuss the global, instantaneous 
character of the dragging of inertial frames as it follows 
from our results and show how it can be seen in diverse 
situations, including cosmological perturbations. We also 
indicate some open questions. Some lengthy expressions 
are written down in Appendices. There, also various inte- 
grals over products of Legendre functions and non-trivial 
radial integrals are evaluated. 



II. ROTATING SCALAR WAVES 

In linearized Einstein theory the basic Regge- Wheeler 
equation [131 ] for odd parity waves becomes the usual 
flat-space 1-dimensional wave equation for the radial 
part. It is thus natural to first present rotating waves 
in Minkowski space. 

Consider the wave equation 



d 2 i> 



■Aip = 



(2.1) 



in spherical coordinates. Assume 
1 



Ira 



[Qi(t,r)Y lm (0, V )+Q l (t,r)Y lm (9,ip)] , (2.2) 



Y lm = N? 1 P t m (cos e)e im *, Z>0, \m\<l; 



(2.3) 



bar denotes complex conjugation, Y\ m are scalar spherical 
harmonics satisfying the standard normalization and or- 
thogonality conditions over a unit sphere (see, e.g., [iij). 



! 2l + l(l-m)\ 
47r (l + m)\ 



(2.4) 



The angular part of Eq. (|2.ip then separates and the ra- 
dial part satisfies 



Qi + Q'i 



2 -q; 

r 



1(1 + 1) 



Qi = Q. 



(2.5) 



Hereafter the dot denotes d t , the prime d r . 

Assume that Qi(t,r) — Q;(w,r)e lwt and substitute 



)i(u,r) = -=ui(u,r). 

\ T 



(2.6) 



Ji is a Bessel function of the first kind. Inspired by the ex- 
ample of the Bonnor [l5[ , Weber- Wheeler [l6| cylindrical 
pulse and by our recent work Q , , in which time sym- 
metric incoming and outgoing rotating waves are smooth 
and finite everywhere at all time, we consider, in spheri- 




Equation (|2.5p is then transformed into 



1 



(l + l) 



h-2 



■ui +uj ui = 0, 



which is the Bessel equation with index v = I + \- 
The independent solutions of (|2.7[) are given in terms of 
spherical Bessel and Hankel functions ji(ojr),ni(u>r) and 
hi(uir). We want solutions regular at the origin, so we 
construct our radial waves as superpositions of 



FIG. 1: Snapshots of function ipi m for / = 27, m = 5, and 
times t — —6, —4, —2, (upper image) or t = 0, 2, 4, 6 (lower 
\"' ') image). For high values of I the wave is concentrated near a 
shell with radius r 2 — a 2 + t 2 , so we can combine waves at sev- 
eral moments into one figure. The wave rotates anticlockwise 
around the z axis. A careful observation reveals that the wave 
comes inwards in the form of a leading spiral, i.e., with the 
outside ahead of the inside (see "white arms" at t — — 2 and 
—4); at t — the spiral structure has changed to a cartwheel 
but rotation keeps the wave away from the origin; at t > 
the spiral becomes trailing (see "white arms" at t = 2 and 4). 
(2.8) At large t the spiral becomes tighter (cf. snapshots at t= 6). 
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cal polar coordinates, the superposition 



,3 r°° 



2r 



(t,r) = Bi\j ^f- / (aoj) l+ 2 e - a " e^JUur) dco 



i+i 



(2.9) 



where the real amplitude B\ and the characteristic width 
of the pulse a are constant. It is convenient to introduce 



a = a{t) = a + it. 



(2.10) 



Employing then formula (8.6.5) in [l7| (see also for- 
mula 6.623 in [3]), we find the integral in (|2 .9[) to yield 



Q l (t,r)=B l 2Hl [(a2+ % 2pi . 



(2.11) 



To construct rotating waves coming in and then go- 
ing out in a time-symmetric manner we just multiply by 
Yi m (9, ip) containing the factor e imv . Clearly, the field is 
the superposition of rotating waves: it is apparent that 
for each lo the wave contains the factor 



(Hereafter, the Greek indices A, /i, v = 0, 1, 2, 3, the Latin 
indices i,j,k = 1, 2, 3; l,m are reserved for harmonic in- 
dices associated with Yi m .) The energy density, the time- 
time component of the stress-energy-momentum tensor 



ST"" 



SC 



reads 



£s = -(^ 2 +7 jfc V^V fc V), 



(3-3) 



(3.4) 



Vfe^/i are covariant derivatives associated with 7^; for 
scalar tp, V kip — dkip, so in spherical coordinates 



1 



</> 2 +^' 2 + 4(cwo 2 + 



r*£ Gin H 



(3.5) 



The energy density satisfies the conservation law which 
can be written as 



(2.12) 



so that the wave pattern rotates with the rate tt p — u/m. 
The rotation and the pulse character of the wave as well 
as its regularity is easily seen when ipi m is written explic- 
itly in real terms: 

B l 2 l l\N] n PJ n (cos9) f l , w N1 
= [((1 + f2 _fi )2+4 fi ]( m)/ 2 cos \™P ~ r ^ < 

(2.13) 

where 



r - t 
r = - , t=- 
a a 



and 



\(t, r) = (I + 1) arctan 



2t 



arg 



1 + f 2 - t 2 
[r~ 2 + (l + ^) 2 ] (i+1) . 



(2.14) 



(2.15) 



Some snapshots of the rotating waves profiles are illus- 
trated in Figure [TJ 



III. ENERGY AND ANGULAR MOMENTUM 
IN ROTATING SCALAR WAVES 



= / V {ipViip)dV = / ipVnpn l dS. 
Jv Js 



(3.6) 



The momentum density pk — —ip'Vkip of the scalar field 
obeys a conservation law as well: 



i I P kdv = i I -^ k ipdv 

at Jv at J v 



(3.7) 



= - / V, 



{TipVkip - 14 [(v^)(v>) - ^ 2 ] } 



\dV 



ThiidS. 



Similarly, the angular momentum in the field obeys the 
conservation law 



-J v e^ P k dV 



e ljk x j T*n s dS. (3.8) 



The angular momentum flux is given by the components 
£ijkX : 'T ks through the surface S. The angular momen- 
tum density along the z-axis is 



We start with some simple properties of scalar fields. 
The energy density of the scalar field can be derived from 
the scalar field Lagrangian density 



(3.1) 



where we use the background Minkowski metric in spher- 
ical coordinates 

9uu = {1, -1m} = U, -1, -r 2 , -r 2 sin 2 9}. (3.2) 



js = -ipd v ip. 



(3.9) 



Let us calculate the energy density and the angular mo- 
mentum density for rotating scalar waves. Consider a 
single (l,m) component, like in (|2.2[) . written now in the 
form 

= \nT [Qi m {t,r)e im * +Q lm (t,r)e- im *] P^) , 

(3.10) 
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with fi = cos 9, Pj m (//) are real Legendre functions. The 
energy density (|3.5j) of the scalar field becomes 

S3 = |(AT) 2 j (QtmQim + QLQ'lm) (^T) 2 (3-11) 
QlmQln 



1-P 



ipn 2 + (i-v 2 )(pry 



i 



(QL + Q' 2 i m y imv + (Q' hn + Q't m )e- 2im<p (PD 



-jm\2 



m 



Qim e ~ 

2 



■(irr + (i-M a )(jr t 



where we write P™' for (d/df^Pf 1 . If we average over ip 
terms, we obtain 



1 

2^ 



(3.12) 



1 



(Nn 2 {(Qi m Q lm + QLQ'i m )(pn 2 



4 

QlmQlm 



1-K 



■(ir) 3 +(i-M a )(jr') s 



If, in addition, we average over the latitude by taking 
j}f-i( € s)d/J., and use the integrals of P"s evaluated in 
Appendix A, we arrive at the energy density, averaged in 
every direction, as a function of r in the form 



16tt 



QlmQlm + QlmQlm + 1(1 + 1 



QlmQln 



(3.13) 



The z-component of the angular momentum density is 
still simpler: 



Js 



— (^D 2 (QlmQlm - QlmQlm^ 



2imip _ 



c.c. )(prr, 

(3.14) 



"c.c." denotes "complex conjugation" . 

So, after averaging over tp and over 9 (see Appendix 
A) we get 



«is» 



1JTL / • * 

Y^Tj; [QlmQlm ~ QlmQlm 



(3.15) 



IV. ROTATING LINEARIZED 
GRAVITATIONAL WAVES 

We consider first the linearized theory of gravity in 
a flat background in spherical coordinates (t,r,8,(p). It 
is natural to decompose the metric perturbations into 
tensor harmonics. These are given in terms of standard 
spherical harmonics (I2.3[) and their derivatives, see, e.g., 



[13j , [l9j , [20j . For any symmetric second-rank covariant 
tensor S^ LU we define the expansion 



Sfi V - ^2 

1, 771 



(4.1) 



+ Ai m (t,r) + Bolm(t> r)boim iiv + Bim(t, r)bi 

+ Qoim(i^)coirn \jlv + Qlm(^-'> r ) 

+ T>im(t,r)di 

771 \JLV 

+ Qim{t,r)gim fj,u H~ J~im(jji ^*)flm fxv i 

where a oim ^(6>, 93), fi miu ,(9,<p) are tensor harmon- 
ics: we follow the notation and normalizations used by 
j20| . The harmonics can be classified into even- and 
odd-parities - those with even-parity transform by the 
parity (—1)' under the inversion, i.e., the transformation 
(9, <p) —> (ir — 9, 7r + ip), in the same way as Yi m (0, <p) 
- while odd-parity harmonics are defined by the parity 
(— 1)' +1 . The harmonics are mutually orthogonal and 
normalized over the sphere of radius r in the following 
sense: 



2-7T f 7T 



'0 Jo 
where 



X Vm ,^P5»°X lmpa sm8d8dp = S n S m , m , (4.2) 



5™ = diag{l, 1, 1/r 2 , l/(r sintf) 2 }, 



(4.3) 



X stands for any of the harmonics a, The orthogo- 
nality properties enable us to find the "coefficients" Aoim , 
etc. in the decomposition (|4.ip of the general covariant 
tensor S^u- For example, in the following section we shall 
need the coefficient Qoi m (t,r) multiplying c ; m . Using 
the orthogonality property of the tensor harmonics and 
multiplying by c / m , we find 



27T p7t 



Qoim(t,r) 



5^6" a S^c i mpa sm9d9di P . (4.4) 



Jo 



It is well known that the resulting field equations for 
the radial factors may be discussed separately for odd- 
and even-parity perturbations of the Schwarzschild black 
hole, see, for instance, 13 1, [13], [HJ, [13], which of course 
remains true for flat space. For our purpose it is sufficient 
to consider odd-parity waves. In the expansions (|4.1[) the 
odd-parity harmonics are Coi m , c/ TO and c?; m , with explicit 
forms of the non-zero components being 



1 



COlm te 



CQlm tip 



C-lm r8 



C-lm reft 



ip *■ Im i 



y/2l(l + 1) sin( 
/ n J, , siu8d e Yim, 

ir 1 n v 

\OpYi 



y/2l(l + 1) sin 
ir 



if * im > 



sin 9df>Y] r 



di 



m ee 



j2l(l + \)(l-\)(l + 2) 



smt 



/ m • 



(4.5) 
(4.6) 
(4.7) 
(4.8) 
(4.9) 
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Mm ifip 



Mm 6ip 



V'2Z(Z + l)(Z-l)(Z + 2) 
V2/(Z + !)(/-!)(/ + 2) 



sin 



sin0W; 



where 



(4.10) 
(4.11) 

(4.12) 



Xi m = 2d v (8g - COt 6) Yl m , 

W lm = (d 2 - cot 9d e - -A^dA Y lm . 
\ sin J 

For the ex plic it forms of the other tensor harmonics, see 
[HL EH; |H| and As mentioned above, our con- 

vention for harmonics is that used in [2(| which differs 
slightly from [l3[ or [lj|. However, following [2(| we as- 
sume the expansion of the odd-parity metric perturba- 
tions h^l ((i) — (1), (2) denotes the first- and second- 
order perturbations) to have the form 



h y m ' m {t,r)c olmfn . 



1 „. li llm\ b i ' ) v lmnv 



(4.13) 



+ i 



V2Z(Z + !)(/-!)(/ + 2) (0 

2 r 2 ll 2lm\ i -' ' ) a lmnv 



so that the factors multiplying the harmonics (j4.5H4.llj) 
cancel the multiplicative factors at h^l, hu , h!$ m in 
(|4.13|) . Hence the radial functions, the h's, are identi- 
cal in [H, pj|, [13 and Q. 

In the literature time-dependent perturbations are usu- 
ally Fourier analyzed and perturbations with even parity 
discussed. In [2l[ a complete account of the Hamiltonian, 
Regge-Wheeler-Zerilli and Newman-Penrose formalism 
for perturbations (with both parities) of a Reissner- 
Nordstrom black hole is given and the dipole I — 1 per- 
turbations are also considered. Putting M = e = 
there and choosing the Regge- Wheeler gauge condition 



(0 



0, we obtain the 



for odd-parity perturbations, h. 

radial functions v (t, r) and v (t, r) in terms of a 
function R satisfying the wave equation 1 as follows (dis- 
carding labels): 



ho = - 



hi=- 



1 



(i-m+2) 

r 

(Z-l)(Z + 2) 



R, 



R + R"- l l±^R = 0. 



(4.14) 
(4.15) 

(4.16) 



We see that last equation can be put in the form of Eq. 
(|2.5[) for the scalar field amplitude Qi(t,r) through the 
substitution 



R = rq. 



(4.17) 



So in terms of q(t, r) the odd-parity metric perturbations 
are given by 



h 



h x = - 



(i-m + 2) 
i 

(*-l)(Z + 2) 



r 2 q. 



(4.18) 
(4.19) 



These formulas are valid for I > 2; in the linear theory in 
vacuum the dipole odd-parity perturbations (I — 1) can 
be transformed away by a simple gauge transformation 
[l9l l2lj |. Later, in Section VI, we shall substitute for 
q(t, r) the solution given in (|2.9[) and consider (|4.18H4.19|) 
to be the first-order perturbations hP-' generating the 
second-order perturbations. 



V. ENERGY IN GRAVITATIONAL WAVES 

In [22j we defined gravitational energy on x°=const 
slices in stationary asymptotically flat spacetimes by tak- 
ing the difference between the total mass-energy and the 
total "mechanical" energy arisin g fr om the matter stress 
tensor. Many formulas given in [22j are true in general, 
non-stationary spacetimes. 

Consider the slice x°=const and the metric in the 1 + 3 
decomposition associated with the Hamiltonian formula- 
tion (see, e.g., [23]) 

ds 2 = goadt 2 + 2g kdtdx k + g ik dx l dx k (5.1) 
= (w 2 - W 2 )dt 2 + 2W k dtdx k - llk dx l dx k , 



i,k= 1,2JL iu = V.900 + W 2 and W 2 = j lk W l W k (cf. 
(2.36) in [22|). The energy of matter can be calculated 
using the Einstein constraint equation for 

2wG m =TZ + K 2 ~ K ik K lk , (5.2) 

where 1Z is the scalar curvature of the 3-metric 7^ and 
the external curvature of x°=const slices is 



(V (i W fc) - |9j7«) . K = <y ik K ik , (5.3) 



Vfe is a 7-covariant derivative. The total (gravitational 
and matter) mass-energy as measured at infinity is M. 
The gravitational field energy on the a;°=const slice can 
be defined as 



1 One can get these equations from the original paper [l3l | as well 
as after going back from their Fourier decompositions in e~ lkt 
to the time-dependence and after correcting several misprints in 
their equations (22). 



E n = M- 



1 

2k 



(Tl + K 2 - K lk K lk )^d 3 x , (5.4) 



k is Einstein's gravitational constant, 7 = dot 7^., all 
volume integrals here and below are taken over x°=const. 
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It is often advantageous to use more general spatial 
coordinates than Cartesian coordinates at infinity, so as- 
sume that the metric of the background Minkowski space- 
time which at infinity is "the limiting spacetime" to the 
physical spacetime, has the general form 2 : 



ds z = dt 2 — r ) ik dx l dx 



(5.5) 



Denoting 7-covariant derivatives by Vfc and introducing 
tensorial quantities 



A' 



ik 



j^m tV"- 
1 ik ~ 1 ik 

^r mn (V l7nfe + V klrn - V nltk ) 



(5.6) 



one can write the scalar 3-curvature density entering (15 .4[) 
as follows: 



(5.7) 



Here C is the Einstein Lagrangian 

C = -Y k (A™ fc A" mn - A m u A n mk ) > (5-8) 
and k % are given by (cf. (3.16) in |23|): 

k % =7~ 1 v fc (ry w ). 



(5.9) 



In an asymptotically flat spacetime the integral of the 
second term on the right-hand side of (|5.7j) gives again 
the mass, 



1 

2k 



d t (^k l )d 3 x 



1 

2k 



^k l dSi = M. (5.10) 



This is the total mass of spacetime as measured at spa- 
tial infinity. In our case of pure gravitational waves the 
surface integral gives the total mass-energy of the waves. 
In a general case of radiation going through spatial infin- 
ity for all times, time-dependent radiative terms ~ 1/r 
would also contribute to the surface integral (|5.10p , how- 
ever for pulse-type radiation considered in this work only 
the static mass-energy results. Returning to Eq in (|5.4p . 
substituting from (15.71) and (|5.10[) we find the gravita- 
tional energy as a volume integral over the slice x°=const: 



Ea = kJ ( KlkRlk ~ K2 ~£)V^ d3 



(5.11) 



with C and Ki k given in (|5.8j) and f|5 . 3[) . The energy 
density sg is a scalar density on x =const due to the use 
of the background metric . 



2 Details on how suitable background metrics and mappings be- 
tween backgrounds and physical spacetime may be introduced 
can be found in [2!. [53 . 

3 It can be shown that (15 ■ 1 1 1 ) is gauge invariant provided that the 
rules of the Brill-Hartle averaging employed in the classical work 
of Isaacson l26t are used (for example, under integrals diver- 
gences become small); in the TT-gauge {eg) becomes ~ hi k h lk 
of Isaacson. 



Assume now the physical metric 7^ to be close to the 
background metric j ik and goo = 1 + /ioo so that 



lik = lik + hik , W k = h k , w = l + -h 00 . (5.12) 

The terms proportional to 0(h 2 ) will be neglected. The 
extrinsic curvature (15.31) then becomes 



Kik = V(ihk)o — „hik 



the Lagrangian (|5.9j) is 



-C = ^7 m hV k h km - \V m hV"h 

- ^\ km v m h lk + \Vh km Vih km . 



(5.13) 



(5.14) 



Since we are considering waves by using tensor spher- 
ical harmonics it is natural to choose the background 
metric 7ifc in spherical coordinates: 



7n 



1 



722 



733 = r2 sin2 ' 



(5.15) 



We start from the metric for odd-parity perturbations 
in the general form (|4.13[) with harmonics given by (|4.5t - 
I4.12[) . However, since in the odd-parity case only four 
(non-diagonal) components are non- vanishing, the fol- 
lowing simple notation will turn out to be useful: 



ho2 — h Q f 2 , ho3 — h f 3 , h\2 — fii/2 , h 13 — h±fs 

--r-zdmYim , / 3 = sin 9d$Yi m . 
smt) 



h 



(5.16) 



Starting from these expressions, calculating Ki k and C 
from (15.13[) . f|5 . 14[) using Vfc in flat-space spherical co- 
ordinates, we arrive, after straightforward - though not 
short - calculations at the following expression for the 
gravitational energy density of odd-parity waves 



e G = ^(K ik K ik -K 2 -C), 



(5.17) 



where 



K ik K it 



K 2 = 



1,2 



2r 4 sin 2 6 



{d v f 2 + deh-2cot6hY 



Ad e f 2 (sin 6 cos 0f 2 + d v f 3 ) 



(5.18) 



and 



-£ = 



2r 4 



ih 1 (rh' 1 -h 1 )(fi + 



+ TT7 ~ deh) 



1-2 * f2 

sin 

2 



(5.19) 



We now express the radial parts of perturbations from 
(I4.18M4.1§|) in terms of the gauge invariant quantity q sat- 
isfying the wave equation (|2.5p (Q there is denoted by q 
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in the gravitational case), and functions /2,/3 and their 
derivatives in terms of spherical harmonics according to 
([5~T6]) . The individual terms in (pU8| to (|5T9| are pro- 
cessed as follows. Each term is quadratic in "hf" , say 
^0/2^1/3- To get real results, we should take 4 

Re(h f 2 )Re(h 1 f 3 ) = ±(h f 2 +^0/2) 5(^1/3 + hih) 
= |(/ioMi/ 3 + hoJ 2 hih) + terms in e ±2m ^. (5.20) 

The last terms drop out after averaging over p - see the 
scalar-field case, equation p.llj) - so we retain only the 
first terms with (j) expressed in terms of qi m (t,r) and 
P ; m (/x). In this way we obtain, after somewhat lengthy 
but straightforward calculations and arrangements, the 
c/j-averaged energy density of odd-parity rotating gravi- 
tational waves: 



2 K (e G ) = l [l(l+l)] 2 Uqij+ q 'q')P 2 



-— 



1-VL- 



■P A + (1 - n z )P 



2r 3 



(r 2 qq)'P 2 



2\ p/2 



+ 



q'q + -m' 



Am 1 



:P 2 + (l-^)P 



m 2 ( m 2 + /i 2 



1 - (j? \ 1 - /i 



1(1 + 1) )P 



--1(1 + 1) ) 11PP' + {p 2 +m 2 )P 
1- M / 




1 - Li 2 



-l(l + l))nPP' 



+ [l+TO 2 -(l-L i 2 )/a + l)l P' 2 



Here P and P' denote Pf 1 ^) and g^P™^); q stands for 

This result has some noteworthy properties. The first 
terms in the first curly brackets coincide precisely with 
the energy (|3.13[) in scalar waves averaged over p. To- 
gether with the next term in [1(1 + l)] 2 they dominate 
for big Nevertheless, for rotating waves when m ~ I 
the scalar field part dominates. The very first term 
~ [1(1 + l)] 2 qq will dominate for high frequency waves 
(qq ~ uj 2 ) and high I for both scalar and gravitational 
waves. 



4 Notice that the terms cannot be commuted since, for example, 
Re(ho f'z) Re(hi fg) ^ Re(ho f^) Re(h\ {2) so that when evaluating 
115.181 1. H5.19H the corresponding terms must be grouped together. 



VI. SECOND-ORDER ODD PARITY DIPOLE 
PERTURBATIONS 

In general the second-order metric perturbations 
can be obtained by solving the equations 

G®[h<»] = -G®[hV,h<% (6.1) 

where represent quantities given in (|4.13p with all 
indices seen explicitly. The right-hand side of (|6.ip is the 
source term in the form of an effective energy-momentum 
tensor quadratic in the first-order perturbations h^> and 
their derivatives known from the solutions of G^J [h^] = 
0. In principle writing down the right-hand side of Eq. 
(|6.ip and solving it is a clear task, in practise a formidable 
one. 

However, to determine the influence of gravitational 
waves on the rotation of local inertial frames at the axis of 
symmetry we do not need to solve (|6.ip in general. First, 
following our method in the translational symmetric case 
@, we shall average over (p so that the source term 
on the right-hand side will become axisymmetric. Before 
averaging it will contain terms of the form dth^ d v h^ 
each of which involves two terms independent of p and 
two terms varying like sin Imp or cos Imp. Those terms 
cannot cause any rotation of the inertial frames on the 
axis. Hence, we shall concentrate on the terms that are 
independent of p and denote them again, as in previous 
sections, by the average symbol (). In addition, as in the 
translational symmetric case the dragging will be fully 
determined by solutions of just one equation in (16. ip 
that with [i = t^v = p. 

Therefore, we consider the equation 

(% p ) [hW] = -(G$[hU, fcW]), (6.2) 
or equivalently, 

R%[h^}=-(R%[h^,hW}), (6.3) 
which follows from (|6.2p because in the background is 

(5.21) 

the flat metric in spherical coordinates and R$ [h^] = 0. 
To solve these equations we shall expand both sides in 
tensor spherical harmonics. Any symmetric second-rank 
covariant tensor can be expanded as in (|4.ip . We are 
interested in the axisymmetric (t , p) component, so the 
only harmonic entering the left-hand side of (|6 . 3[) lin- 
ear in is the odd-parity harmonic coi m given in 
(|4.6p . In order to extract the contribution proportional 
to the same harmonic on the right-hand side of (|6.3p we 
imagine it to be expanded as the general covariant ten- 
sor in (|4.ip . The coefficient Qoim(t, r ) multiplying 
coim^u can be derived using formula (|4.4|) . In our case of 
the averaged - and thus axially symmetric - component 
(R < j. 2 J[h^ 1 \ h^]) we put m — and (fiv) — (tp) to obtain 



Qow(t,r) 



2n />7r 



(2), sin 9d8 

-2{RL ) c oiot v 9 . 2 „ dp 





47T 



r 2 sin 2 1 



y/2l(l + l)r Jo 



(R$)d e Y l0 d9, (6.4) 
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where we substituted for coiotip from (14.61) . Hence, the 
right-hand side of equation (|6.3|) reads 



2 J {R^[h^,h^})d e Y lo d9 



sin fldgYjo 



(6.5) 



The left-hand side of (|6.1[) in flat background (in spher- 
ical coordinates) is formed from the second derivatives of 
the second-order perturbations: 



h (2);<x 



(2) ;a/3 
/3a 



A(2)ia _ 
l(2)0;c 



(6.6) 



We expand G^J in tensor spherical harmonics and realize 
that we need only the axisymmetric component ; the 
result will thus be proportional to co;o- It is now impor- 
tant that by using the gauge transformation x' — > x+£ to 
the second order we can choose the Regge- Wheeler gauge 
also for the second-order perturbations (HI, so for 



(2) 



and 



odd-parity perturbations we can again set h 

keep just and 

Calculations of G& [h^] from (I6.6[) lead us to the fol- 
lowing expression: 



. /l (2)" + /G + l) /i ( 2)+A ( 2 ),. 



2- 



(2) 



sin0d e Y w . (6.7) 



Combining this with the right-hand side (|6.5p we arrive 
at the differential equation implied by the field equation 
(|6.2|) for radial parts in the form 



l{l + l) u {2) ;(2)/ 2,(2) 

I2 ,L o 'n ~ ri i 



(6.8) 



2tt 



1(1 + 1) Jo 



(R$[hW,hW])d0Y m dB. 



This equation suggests that we need to consider more 
field equations to determine both second-order pertur- 



bations and h\ z> from the first-order ones. In- 

(2) (2) 

deed, if we wish to calculate both and h\ ' for all 
Vs we have to take into account equations (|6.1[) also 
for (fJ/V) = (9, if), {r,tp) and by combining them we can 
achieve just one wave-type equation for h( 



,(2) 



(2). 



8C 



■,(2) 



— S ri p + S'g v Se<p, (6.9) 



where components S^ v stand for the averaged compo- 
nents {R$[hW,hW]). This is a hyperbolic equation in- 
dicating that for general / the effects of the first-order 
terms are perturbations of the second-order terms that 
are non-instantaneous. However, in all our previous work 
on dragging of inertial frames due to angular momen- 
tum of the sources (see, for example, Q and references 
therein), the effects were instantaneous, at least at this 



lowest order. In fact, the same situation arises here. In- 
ertial frames at the origin will be influenced primarily 
by the dipole perturbations since all waves behave as r l 
there. Now it is well known that for I — 1 one can achieve 
hi = by an appropriate gauge transformation [l3l |. [2lj |. 
For the dipole second-order perturbations the equation 
(16.81) becomes 



,(2)" 



-h (2) - 



(6.10) 



where 



9(t,r) = -d — 




sin6d6dip. (6.11) 



Here we indicated the </j-averaging explicitly (which can- 
celled the factor 2ir) and we substituted for deYio- With 
g(t,r) known, the solution of (|6.10[) can easily be found 
by variation of constants from the simple solutions pro- 
portional to 1/r and r 2 of the homogeneous equation. 
Putting the arbitrary constants entering the homoge- 
neous solution equal to zero, the complete solution for 
h^\t,r) reads 



,(2) 



■± j o 9(ty)r' 2 dr' - T - J 9(ty)^. (6.12) 



VII. HOW ROTATING WAVES DRAG 
INERTIAL FRAMES 

Since we are interested in the dragging of inertial 
frames near the origin we wish to express 



9(ty)±jdr', (7.1) 
o r 



with r) given by (|6TI|) and (R$) in Appendix B. 

The general form of the odd-parity metric perturba- 
tions (|4.13p . with the tensor spherical harmonics (|4.12l) . 
for the special case of I = 1, m = perturbations with 

(2) 

only h 10 ^ yields the metric component g tip near the 
origin in the form 



(2) _ 

9 tip ~ 



'—h 2 \t,r) sin 2 6 

47T 



-uj q r 2 sin 2 9, (7.2) 



where the rotation of an inertial frame (of a gyroscope) 
located near the origin has the angular velocity 



w 



1 

47T 

sin( 



dn 



dr 



(7.3) 



where dfl = sin 8 d9 dip. 

To evaluate [h^\ h^}, which has more compli- 
cated structure than any expression calculated so far, we 
found it more convenient to start from the real form of 
the solution of the original scalar wave equation (|2.1I) and 
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use it to express the gravitational wave perturbations - 
see Eqs. (|7.4j) . (|7.5|) below. As a consequence, also the 
metric coefficients will be real from the start, rather than 
of the form (|5.16[) . With real expressions one does not 
have to be careful when grouping together properly the 
corresponding terms in the quadratic expressions like we 
had to do in calculating the energy densities (cf. foot- 
note 4) and it is much easier to find the results using the 
standard (commutative) computer algebra codes such as 
tensor package in Maple. In fact, we first proceeded by 
using the same methods as in the case of energies and, 
after going to the real variables, we checked that both 
results agree. 

So now we start from the metric perturbations written 
as follows 

1 Qx Qx 

hte = ^9d^> ht * = - S[n9 frd9 > (7 ' 4) 

1 d~v d\ 
Ke = . Qfua , h rv = sin9—-— , (7.5) 
sin 9 otdip otoO 

where function x(^ r i$i¥') = ^ V'imj (ip satisfies the 
scalar wave equation (|2.1j) in flat space); introducing x 
instead of ij) turns out to simplify the second-order Ricci 
tensor. Using the general result for the the component 
(see (|B2|) in Appendix B) and substituting for \ 
from (17.4117. 5[) we arrive at 



average over the spherical angles we first find 



R t 2 J = -| ^2 \^Xfirr8X, V t ~ X,r<prX,9t$ (7.6) 
+ 2 X,tr V 8X,r8 — X,rr<p8X,t6 ~ Xfitr8X,ipr + X,rr9X,tip9 
+ X,Br9X,ript ~ X,t8X,tt<p6 ~ X,tt9X,ttp6 ~ X,^ptX,9US^j 
, {X,rBX,r V t ~ X,tr6X, V r ~ X,r V rX,t6 + X,tt6X,<pt) COS 9 

2r 2 sin 9 

X,r(f(ptX,<pr + X,tptprX,rtpt — X,tiptX,(p<pt X,r8X,<pt COS 9 



1 1 



2r 2 sin 2 

X,6reX,yt . X,ft 



(2x,v»pr + 2rx ,(ptptt TXwprr) • 



r 5 2r 3 sin 9 

We shall now write the real part of the scalar potential 
x(t, r, 9, ip) in the form corresponding to the form (|2.13[) 
as the product 

x = BiN™ K(t, r)P l m (cos9) cos(mip - X(t, r)) , (7.7) 

where Bi = Bi2 l l\, N™ and A are given in (|2.14j) and 
(|2.15|) . and we have introduced function 



«(t,r) = 



[(l + r 2 -t 2 ) 2 + 4P]('+ 1 )/ 2 



(7.8) 



Inserting the function x from (JT77J) into (|7.6I) and em- 
ploying the formulas for the integrals of the products of 
the Legendre functions given in Appendix A to find the 



R[ 2 J dtt - -B?^ 



+2 (l 2 + I + l) KjK^r^.r I (I + 1) K (X t rK>,t ~ ^,rA,t) 

+ (I 2 + I + 2) K (n, t \rr ~ \tK,rr) -1(1 + 1) K 2 A, t 3 

+1 (I + 1) (k 2 A,* A,,. 2 - 2 ft t 2 A :t + KKjtM - K\ttK,t) ■ 

(7.9) 

Substituting for n and A the expressions (|7.8|) and 
(|2.15[) and integrating over the radial coordinate we ob- 
tain the angular velocity ujq in the form 



Bf m(l + l)(l + 2) 



Mt) = tt ,, iqw 
2ir a (t + 3) / 



(7.10) 



Ui-Vif 



l + t 2 )l? 



l+3 (t)+[U l + V l t \I{ +3 (t) 



where 

Ui = (2 ^ 5 + 7 / 4 + 4 Z 3 - 7 / 2 + 24 / + 36) , (7.11) 
Vi = 3 (Z 4 + 2 Z 3 + 3 Z 2 + 8 Z + 12) . 

Here the radial integrals iff are defined and analyzed 
in Appendix C. Although the last result has indeed an 
analytic, explicit character, it is not very practical for 
calculations because the integrals J;_L 3 are determined by 
high-order derivatives for large I. Hence, we describe 
their properties in detail in Appendix C. 

Turning to the numerical evaluation of the angular ve- 
locity of dragging we obtain nice "bell-shaped" figures 
for ujq for various values of m at a given I. In Figure [2] 
ljq is plotted for / = 10 and m = 1, 2, 3, ...10. 

Also, the integrals can be explicitly calculated at the 
moment of time-symmetry when the dragging near the 




FIG. 2: Angular velocity of the central inertial frame 
uJo(l,m;t) for I — 10 and m — 1,2, ...10 (from bottom to 
top). The vertical axis is scaled in units of wo(10, 1;0). 
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FIG. 3: The dependence of normalized angular velocity of the 
central inertial frame u)o(l, l;t)/uJo(l, l;t = 0) on the parame- 
ter I = 2, 3, 10, 20, 30 (from inside to out). Also the function 
(1 + P)~ 3 ^ 2 is shown as a dashed line to indicate the limit for 
large I (cf. Eq. (Eg) ). 



which describes a ray from the star with celestial coor- 
dinates 9, ip which reaches the center at time T. The 
geodesic equation can be written in the form 




- 2r<°>* (x v = x(°)"(A)) pWPpW' . 

To solve these differential equations we use covariant 
components of the perturbed photon momenta for which 
the effects of curvilinear coordinates cancel out; we omit 
superscript W for simplicity. Then we can write 

T\ p ° ) = ~ r2 lx p(1)e - 2rpil)6pr (8 - 4) 

= -(p*) 2 hte,t - {p r ) 2 h r e tr - p t p r [h re ,t + h te ^] 
4tP { ^ = -r 2 sin 2 e4rP {1)v ~ 2rpW* sin 2 9p r (8.5) 
= ~ip t ) 2 h tv>: t - (p r ) 2 h rVir -p t p r [h rtp ,t + ht Vt r] ■ 



origin is maximal since the rotating waves are at nearest 
distance. The integrals have simple forms also at t = ±a, 
but these values are of no special significance though they 
can serve for checking the numerical results. The central 
and maximal value reads 



wo(t = 0) 



Bfm{l + 1)\ (Z + 2)! 
4777, Z (2 Z + 3)! 



Ut 



(7.12) 



Finally, at late times the analytic behavior of I L 
plies 



M 



wo(t = 0) 



a 


•w 


a 
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7 





(7.13) 



and in Figure [3] we illustrate that ui(t) for large I and 
given m can be approximated by function 



[1 + i 2 ] 1 



(7.14) 



VIII. 



OBSERVING STARS THROUGH 
GRAVITATIONAL WAVES 



Since our waves are sufficiently weak near infinity, we 
can consider static stars there and ask how they would 
move across the sky as seen by an observer fixed in the 
flat region at the origin. We shall evaluate the first-order 
effects of the waves on the propagation of photons, the 
effects due to dragging of observer's frame are of sec- 
ond order and will be discussed subsequently. The influ- 
ence of gravitational waves on the passing photons can 
be computed as a perturbation of the ingoing radial null 
geodesic 



,(o)m _ 



= [T + A, -A, 6, <p], A € (-oo,0] , 



(8.1) 



We see that along the radial rays (for which sin 
const) the covariant components of photons momenta 

Pg 1 ' and p^ can be written as integrals along the 
unperturbed ray and we do not need to solve the geodesic 
equation as set of differential equations. Similarly, the 
differential equations for 9^ and ip^- 1 ' can be written as 
follows 



d 



d (1) 



1 (i) 



r 

LA 

j/dA 



lid 

p r r dX^ 



.(!) 



1 



_p(i) 



(8.7) 



1 d 

pr~d\ 



r sm 



_„(i) 

2 q^V 



1 



1 



p r r sin 2 9 dX 



__p(i) 
2 a A \ *V 



For photons which reach the origin the boundary terms 
arising after integration by parts vanish, so we can write 



50{T) = 



Y-Tx p( M 

p r r dX 



1 1 d 



o (P r ) 2 

[Ke,t 4 



rdX 

ht9,r 



Pgdr 
hte,t 



(8.8) 



dr 



i\t=T- 



6<p(T) = 



1 1 d (i)^ 
~ ■ 2/, 7TPy dX 
oo P r sm 9 dX v 



(p r ) 2 rsm 2 9dX^ 



(8.9) 



dr 



L rtp,t 



Hip,r 
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Since now we work in the linear approximation to find the 
light deflection we can go back to the simpler, complex- 
valued formulas (|2.11[) and (|4.14M.19|) . When we intro- 
duce amplitudes 



A9 

Aip 



mP™ (cos ( 
sin 8 



-BiNT2{l- 1)!P7(cos0), 
the integrals (|8.8[) and (|8.9|) simplify to 
59{T) = Re[A9 p(a + iT)], 
5ip(T) = Re [Aip ip(a + iT)] , 
where the complex function p(a + iT) reads 



-i e 



,1(1 + 1)0 + 2) 



J+2 



(2r) 



1-1 



(a + iT) 1 1 Jo (a + iT ~2ir) l+3 



(8.10) 
(8.11) 

(8.12) 
(8.13) 

dr 



1 



. T \l+2 



.14) 




If we assume a telescope pointed towards the star's ini- 
tial position at T — > — oo, the quantities 89 and sinOSip 
are proportional to star's coordinates on telescope's pho- 
tographic plate; star's trajectory can be described by a 
simple formula 



Sip . 59 
~Aip~ + l ~A9 



i l e i 



1 



(8.15) 



We see that the same image is replicated for all stars 
(Figure 0]). Depending on the star's position the image 
is rotated by e tmLp and then scaled by factors A9 in the 
latitudinal and by Aip in the longitudinal direction. Sur- 
prisingly, the images are time-symmetric with the max- 
imal deflection occurring at T = despite the fact that 
the deflection is calculated using retarded integrals (|8.8[) 
and ASH) . 

Since the calculation of the second-order metric every- 
where would be very complicated we cannot, of course, 
describe all second-order effects. However we can de- 
scribe the second-order effect of the total rotation of the 
central observer's telescope due to the dragging of his in- 
ertial frames given above, cf. (|7.10l) . (|7.11l) . Since the 
pulse starts and ends at infinity at t = ±oo, the total 
rotation of the telescope is clearly gauge invariant. As 
a consequence, the curves in Figure 0] will no longer be 
closed. The stellar initial position at T — > — oo will dif- 
fer from the final position at T — > +oo by an amount 
proportional to the area under the graph of u>o(t). 



1 = 3 1 = 13 




1 = 23 1 = 33 



FIG. 4: Since light from distant stars is influenced by the 
gravitational waves the observed positions of the stars change. 
An observer at the origin can record the apparent position 
of the stars on the celestial sphere on a photographic plate. 
When appropriately scaled and rotated, the trajectories of all 
stars are the same. A star starts at the origin of the plate 
(x = y = in the planes above) and moves along closed 
trajectories the structure of which becomes more complicated 
with increasing /. At the moment of time symmetry, t = 0, 
the image is located at maximal value of x (with y = 0). 
This curve is an image of a straight line in the complex plane 
z = 1+is, s € ( — oo, oo) mapped by the function f(z) = z~ l ~ 2 
(cf. Eq. I|8.15|l h Together with the trajectory, positions at 
time t = 0, ±0.05, ±0.1, ... are shown as circles. 



IX. CONCLUDING REMARKS 

The results obtained here extend significantly our past 
studies 0, i| of gravitational waves with the translational 
symmetry that rotate about the axis of cylindrical sym- 
metry and cause a rotation of local inertial frames at 
the axis. In the present work we investigated the effects 
of the rotating gravitational waves in a general, asymp- 
totically flat situation. The main conclusions following 
from these studies demonstrate explicitly that any state- 
ment of Mach's principle that attributes all dragging of 
inertial frames solely to the distribution of energy and 
momentum of matter, characterized by T^ v , as the ori- 
gin of inertia is false - gravitomagnetic effects are also 
caused by gravitational waves satisfying Einstein's field 
equations in vacuum for which a local energy and mo- 
mentum cannot in general be defined. 

We found that the dragging of inertial frames at the 
origin occurs instantaneously \ with no time-delay de- 
pending on the position of the pulse. The dipole pertur- 
bations determining the dragging follow from the con- 
straint equation which is elliptic (see equation (|6.10[) h 
Gravity is a tensor field, so gravitational waves cannot 
have a dipolar part. This is the case with our first-order 
rotating pulse. It is the non-linearity of Einstein's field 
equations causing a self-interaction of the waves which 
produces an effective source for the second-order per- 
turbations. For a general multipole moment I > 2, the 
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resulting equations are again wave-type equations indi- 
cating that the effects of the first-order terms cause the 
second-order terms in a non-instantaneous way. How- 
ever, inertial frames at the origin will only be influenced 
by the dipole perturbations since waves decay as r l at the 
origin. And for the dipole (I = 1) second-order pertur- 
bations, the equation becomes elliptical, with the source 
term given by a component of the second-order Ricci ten- 

(2) 

sor R\' quadratic in the first-order metric and its deriva- 
tives. 

The instantaneous character of the dragging of iner- 
tial frames by angular momentum of sources implied by 
the Einstein constraint equation can, in fact, be seen in 
diverse situations, including those in cosmological pertur- 
bation theory. We encountered the effect in the study of 
local inertial frames affected by the rotations beyond the 
cosmological horizon Q, [lj|, or in the general analysis 
of "rotations and accelerations" of local inertial frames 
in the linearly perturbed Friedmann-Robertson- Walker 
spacetimes Q . In the cosmological context we introduced 
desirable gauges - we call them the "Machian gauges" - 
as those in which local inertial frames can be determined 
instantaneously from the distributions of energy and mo- 
mentum in the universe by means of the perturbed Ein- 
stein constraint field equations. The "uniform-Hubblc- 
expansion-gauge" or the "minimal-shear hypersurface" 
gauge combined with "the minimal distortion" gauge on 
the spatial metric are examples of the Machian gauges; 
the synchronous gauge or the generalized Lorentz-de 
Donder (harmonic) gauge are not. (See [|[ for detailed 
analysis and preferences of the Machian gauges.) 

Likewise, in the context of gravitational waves we 
should bear in mind that such physical effects like the 
dragging of inertial frames are of a global nature and re- 
quire the introduction of suitable coordinates in which 
the waves are described 5 . 

In the present work we also investigated for the first 
time, as far as we are aware, the effects of the rotating 
gravitational waves on the propagation of light from dis- 
tant fixed stars at (asymptotically flat) infinity. 

There are various open problems in these themes which 
deserve further study. In order to have more natural un- 
derstanding of the effect of rotating gravitational waves 
one should consider the waves within a cosmological (not 
asymptotically flat) context. Although in [3j we empha- 
sized the role of the instantaneous Machian gauges, it 
would be also interesting to study these effects, in par- 
ticular those connected with gravitational waves, in other 
type of gauges like the Lorenz-de Donder gauge (cf. Ap- 
pendix D in Q). 

Of course, any simple model problem illustrating the 
Machian effects exactly within general relativity, like that 



5 This point has been emphasized also in contributions and dis- 
cussions in the excellent book Q (see, for example, remarks by 
D. Brill, p. 213). 



employing the conformal stationary metrics in our recent 
note on the linear gravitational dragging [27]], deepens 
considerably our understanding of such effects. 
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Appendix A: Integrals over products of Legendre 
functions 



The following integrals over 8 or fi = cos 8 are needed 
to evaluate functions K(t, r) in (|6.12p in terms of radial 
integrals over g(t,r) given by (|6.11[) and (|B3p . Most of 
them can be found in |18j . two follow from simple inte- 
gration by part. These are (0 < m < I): 



x 2/ + 1 (t — to)! 

_i 1 — /i 2 m (I — to)! ' 



(i-^)(JT')> 



(l + m)\ 
{I - to)! 



21(1 + 1) 
21 + 1 



(Al) 
(A2) 
(A3) 

(A4) 



Appendix B: Calculating and averaging the 
second-order Ricci tensor component 



To compute R^[h^,h^] we use metric perturbations 
(|7.4p . (I7.5P which in coordinates — {t, r, 8, ip} yield a 
linearized metric with odd-parity perturbations on a flat 
background. We start from the general form 



Qia> 



1 



hte 

htip 



hte 

-1 h r e 

Ke -r 2 

Ke 



Http 

h<r(p 



- 2 sin 2 8 



(Bl) 



With the metric functions given by (|7.4I) , (|7.5[) the com- 
ponents of the Ricci tensor in the first-order approxima- 
tion yield zero either identically or due to the fact that 
ipim solves the wave equation (|2.ip . The only second- 
order component we need is 



— R 



(2) 



ht v fiK 

fh r g g 



htip : rhr 



+ Kg(h tip jg + h te ^ t ) + h r ip,t' 

+ Ke(hte,ipr + h rip je — 2ht<p, r e ~ 2h r g iVt ) + hte,tph r 6,rj 
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Effects of rotating gravitational waves 



1 cos 

— t: (httphtQ t — h rtp h t 0. r — htehtui t) 
2r z sin 



(B2) 



1 1 cos 6 

H — nh r e (ht v ,e — h t e, v ) + 7r^~ — nhr8 (h rv ,t + h tip , r ) ■ 
r d 2r z sm8 

This simplifies to (|7.6[) when perturbations in the form 
(|7.4|) . (|7.5|) are introduced. When integration over the 
spherical angles ip and 8 is performed, we get (I7.9[) . At 
this moment we use functions k and A given in (17.81) and 
(|2.15p and obtain 



. dr 



2B?m(l + l)a 2l+5 r 21 - 1 dr 



4a 2 t 2 



1+3 



| r\ (ri + 4 a 2 t 2 ) (l 4 + 2 / 3 + 3 1 2 + 12 1) (B3) 
+ 4 (t 2 + a 2 ) [3r j - 4 r 2 (2 a 2 + 2 r 2 + t 2 )] 
+ 8 r 2 Z [a 4 + r 4 + 2a 2 < 2 - t 4 - 2 (/ + 4) a 2 ^] j , 

where r 2 ^ = r 2 + a 2 ±t 2 . After the numerator is expanded 
in powers of r, using t = i/a, this integral can be con- 
verted into a linear combination of integrals Iff defined 
in Appendix C: 



Za? +3 + (Z, +8l + 12) (r 2 + l) 3 7f +3 (B4) 
+ (t 2 + 1) [3Z, - 16(7 + l) 2 + 8 - (Zj - 8) t 2 ] /, 2 +3 
+ [3Z ; - (41 + 5) 2 + 5 - (Zi + 8 / + 20) r 2 ] // +3 , 



where Zi = I (I 3 + 2 1 2 + 3 1 + 4) . Using identities (jCl~5)) 
in Appendix C the number of integrals iff can be reduced 
down to two, as is shown in the final form of the formula 
for uj in (|7TTU|) . 



Appendix C: Radial integrals iff , jff 



To evaluate the metric function ho at a general radius 
r we need indefinite integrals in (16. 12[) ; to determine the 
rotation of the inertial frame at the origin we have to 
evaluate the definite integral in (|B3|) . If we introduce 
dimensionless variables £ = r 2 /a 2 , and t = t/a, the defi- 
nite integrals we need to evaluate simplify into 



c 



L-M-l 



T d C (CI) 



h [(l+T 2 ) 2 +2(l-T 2 )C + C 2 j 

where L will be I — 3, M is a small integer. We first put 



c 

1+T 2 



P 



1-T 2 
1+T 2 



(C2) 



so that 



(1+T 2 ) 



2\L-M /-oo 



„L-M-1 



(l + r 2 ) 2L J (l + 2fix + x 2 ) L 



dx. 



(C3) 



To evaluate this integral we start from the following sim- 
ple integral (with a=const parameter) 



H 

Setting then 
Hm = 



dx 



o a + 2(3x + x 2 



(C4) 



arctan 



dx 



a + 2j3x + x 2 ) 



2\M+1 



d \ M H Q 



da / M\ 



, (C5) 



we can find the integral which appears in (IC3[) by evalu- 
ating (L — M — l)-th partial derivative of Hj\f(a, (3) with 
respect to /3, so the original integral turns out to be 



t m _ 

It — 



2 M (1 



J2\-L-M 



d L - x H a 
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Since it is inconvenient to use this general form we will 
write the two integrals which appear in (|7.10[) explicitly: 



rii 
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(r 2 + l) 3 



arctan r 



(r 2 + l) 2 d 
t dr 
- 4 -l 
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4 -l 
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(3r 
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(C7) 

(C8) 
3) • 



The expressions for these integrals become much sim- 
pler at t = (and at r = ±1 which, however, we will not 
discuss further). At the moment of time symmetry the 
denominator of the integral (jCll) simplifies to (1 + £) 2 , 
and we get integrals of the form 



rM _ 
j l — 



(i + xy 



-dx. 



(C9) 



Assuming L — M > 2 (which is satisfied in the cases we 
need), we can easily integrate by parts and get 



,M 



L- 1(1 + SB) 



L-l 



M 
L-l 



r M-l 

'L-l ! 



(CIO) 



where the first term vanishes if L — M > 2. Continuing 
with integration by parts we arrive at the simple result: 



7 M 



M\(L — M — 2) 
(£-1)! 



(Cll) 
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This yields Even though (|B4|) contains integrals lM 3 with M = 

0,1,2,3, only two of them are "independent" because 
tM ( = ()) = (£-M-l)!(L + M-l)! of the fact that J°° = for functions satisfying 

Ll ; (2L-1)! 1 ; /(0) = /(oo) implies 



For large |r| dC2j) and ((G8J provide 

(2Z + 3)!! 7T 
2 3 <+ 9 (/ + 2)! |r| 3 



'/ - - Si + O (05) • (CIS) (L - M - 1) (1 + r 2 ) 2 If + 1 (C15) 

- 2M (1 - r 2 ) /f - (L + M - 1) if - 1 - . 



'^S^y.W + iw l ' f ' 1Ji 
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